General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 



A Reproduced Copy 

OF 




Reproduced for NASA 
by the 

Scientific and Technical Information Facility 


H. G. ELROD, Jr. ! 

Profei^or, 

Depcr’men? ef Mechanical Engineering, 
Michigan State University, 
East Laming, Mich. 
Assoc. Men-?, ASME 



deg 4 wy* 
Lewis library, nasa 

CLtVElAN,'), L*HI() 


- * ■ P.tip'sr No. 

62 — WA -221 

The Theory of Pulsating Flow in 

Conical Nozzles N 64 81899 

(JWzlX-J /J ° 

.] knowledge of the dynamic charai terislii s of nozzles and orifices is important in many 
coal ml and stability analyses of engineering ilcviirs, It is usual In assume that the in - 
stanlnncoiis flowrate, for a given set of inlet conditions and. outlet pres sure,, is the. same as 
the nontransient value jin the same operating t auditions. Recently, in connection with 
the stability analysis of an externally pressurised thrust bearing, the Validity of this 
assumption was questioned. The analysis presented in this paper was undertaken to 
phndde an answer. The present analysis applies to any /laid, liquid, or gas f owing 
into a simple eonieal nozzle. The amplitude and phase of the mass-dux response to a 
sinusoidally time varying pressure JI initiation j it the nozzle exit are determined. An 
tip pi oximate formula is given for these quantities in terms of the nozzle throat coca, the 
solid angle subtended by the nine, the. velocity of the lluid ot the nozzle throat, the acoustic 
veloi ity at the throat, and the frequent y of the pressure t’lu tuntion. 


Noz/t.r.s and orilitvs have innumerable ii]>jili<>iil.iotis 
in engineering devices. Freq.u nl.lv a knowledge of the pep. 
formatter of such nozzles and entices under dynamic conditions 
becomes important. For example, on.’ may iio concerned with a. 
control or a stability analysis in which wane relation between Mow 
fluctuations and e\it pressure fluctuations must ho incfirporatcd. 
Tin* usual analytical procedmo in those circumstances is to 
postulate a iplnsi-st e.tdy relation ltd ween flow' rate and exit 
pressure: that is, the tlmv rate is taken, at every time instant,, to 
correspond to the instantaneous exit pressure. This procedure 
cannot lie adopted without some misgivings as to its validity, 
especially if the frequency of the exit pressure llnct nut, ions is 
high. 

Because of the importance of flow llnel nations through nozzles 
and orifices, considerable literature on the subject, has de- 
veloped. Much of t his literature has been summarised in a, survey 
article by Opponhoint |lj. 8 The majority of investigations has 
been experimental in nature, and most of these investigations have 
been devoted to determining the effect, of pressure (hiol nations on 
otrroirf flow . However, Sehiill.t-tirunow j ti | repot (s that if the 


1 Also, Consultant, The Franklin Institute, Philadelphia. Pa. 

8 Numbers in brackets designate Itefeieiiees at end ot impel. 

Donfrihuted by the Applied Meeti.iuic's Divi-ioii foi preseitt.at itm at 
the Winter Annual Meeting, New \ orU, N. Y.. N.neather git HO, 
llHtS, of Tup Ami title tx Socikti .ot- M I'.ettA Nte.M. Mnoinuuks. 

. Discussion of this patter should he nd.hess.cd to the Ftlilorial De- 
partment, ASMK, li oiled I'ligiueering ('enter, .'it. r ) Must I7lh St, reel , 
Now York 17, New York, and will It.' accepted until one month after 
final publication of the paper itself in llte JofPiNAi. of Aeei.n, to 
Mf.OHvnh’S, Manuscript received b,\ ASMM Applied Mechanics 
Division, January.), Bh'2. Pspci No Ok# - WA-221. . 


so-rajled “Stroiih.il number," S - ixD/tl, is less than 0.01, qttusi- 
steady pressure-flow relalions are valid. (Here W is the angular 
velocity of the duel tinting exit pressure, I) is the throat diameter, 
ami l r iH the lluid velocity at t he throat,.) Corresponding theo- 
retical work appears lo have been devoted principally to a study 
of the behavior of individual pulses, rather than of frequency 
response. •’ 

The purpose of 1 he present work is lo predict the amplitude and 
phase relations between the instantaneous flowrate and a sinu- 
soidally llnel uating exit, pressure in a conical nozzle, Fig. 1. A1-, 
though the chosen geometry is, of course, extremely simple, com- 
pared with many in use, it is believed that the essential churac- 
terisl it s of a compressible lluid expanding and accelerating from a 

1 A paper using reflection tomb sis lo deal willt problems' akin ,to • 
(host' It riilt'd liete is that of I'.wvcll If) k 



Fig. 1 Conical noxilfl thawing converging ipherically tymmelric 
inflow from stagnation conditions at infinity 


Nomenclature — — — 

o w local acoustic velocity 
.1 »= local (low area 
J * =- flow a rett where .1/ -= 1 
T » amplitude of oscillation 
it <*■ fluid enthalpy 
it m stagnation . m'aalpy. h (•*.'# 

X ~ dimension!) ss perturbed stagnation enthalpy 1 * U'/uy 
M = local Mach number, i fa 
X ~ amplitude factor 

p — fluid pressure 

P •« umplitudt'of p’/p oscillations 

() "■ tunplil, title of ( pv)’/(pe) oscillat ions 
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r - tradius front cone apex 
r = also, residual function 
It = rtidius, r, where M •- 1 in conical flow 
.S’ fluid entropy 
S = also, Ktrouhnl immher, ut)/v 
S'* = mndifled Ftrouhal number, jivl /o 

I - time 

T -■ absolute temperature 
U - fluid velocity ;it nozzle throat 
r ■- local fluid veloi ity ( vectorial in hold face) 

V auxiliary function 

te — dimensionless angular velocity of oscillation 

1 



stagnation reservoir :ire well enough realized for the results to 
have semi-quantitative practical significance. It. is recognized 
til, at near the aperture of any eonie.al nozzle the flow will rouse lo 
be spherically symmetrical but, attain, this effect is not. deemed 
large enough t.o defeat the purpose, of the investigation. 

In the analysis, the flow is presumed to he adiabatic, and frui- 
tion-free. Such assumptions are well known to ho in close accord 
with the facts, f'nder such eireumstanees, Kelvin’s theoreni 
also assures that the flow is irrotational, since it originates in a 
reservoir of stagnant fluid. The fluid itself is assumed to exist, 
in thermodynamic equilibrium and in a single phase. Otherwise, 
the nature of the fluid is praeticalh unrestricted; i.o., it. may bo a 
liquid, vapor or gas. 


Basic Propagation Equation 

To start, the present analysis, wo write down Killer’s equation, 
the mass continuity equation, mid the so-called “TdS" equation 


of thermodynamics. Tims 
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Here D/Dt — d/0 1 + v-V is the Stokes or J.amb operator, and 
the other symbols, which are standard, are defined in I, lie Nomen- 
clature. 

For an irrotational flow: 

(W'v - V '.j ( l) 


so that equation ( 1 ) can he rewritten a 
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But since also, from equation (A) 
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Xh - Vp 
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in an iscut Topic flow, we obtain In mmhining equations (f>) and 
t(i) the following useful result; 



(7) 


where II is the stagnation enthalpy; i.e,, II •-> It + r 2 /2. 

Another result involving the stagnation enthalpy is obtained 
by first, taking the scalar product of equation ( I ) with the velocity 
v. Then 


— Nomenclature— — - — 

O! = ratio of local velocity to acoustic velocity fur upstream 

= ratio of acoustic velocity far upstream to local acoustic 
velocity, «../« 

y, — local fluid isentrujiie expansion exponent, (0 In p/0 In p)s 

y* -■ dimensionless fluid eoeff;eient, 

y * - 1 - (0 '0 In />,)., In (O/i /'dp \ 

f “ ratio of local radius to that where M ■” I in conical flow, 

r/lt 

i? — .’Hixilian function 
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where use lias iiecii made of eipiation (0). Finally, rearrangement 
of equation (Mi, together with the use of the isontropie relation 
between pressme and enthalpy, gives 


1)11 I Op 
1)1 *’ p Of 


((» 


Now the velocity of sound in an isentropie, medium is 

» -. |f0p/0p)s]’'' ! (10) 

As a eolise<|Ueiiee, we may write 

1 !>! I 1 Op 0 In p 

<r- hi " p of ~ ’of ' 1 


Taking the Stokesian derivative of equation ) 1, we obtain 


D 1 DU I ) 0 In p 
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But: 
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Substitutions into equation ( hi) from equations (2) and (7) yield 
the following “propagation equation” for the stagnation en- 
thalpy; thus: 


I) 1 Dll 

lit 1 1 2 ni 


x-u V//-V In p 


(H) 


Differential Equaiion for Spherically Symmetrical Small 
Disturbance 

Let, us now suppose that the flow through the nozzle deviates 
only very slightly (and in a spherically symmetrical manner) 
from a steady, conical inflow. Moreover, in what follows let us 
denote as primed quantities all fluctuations from steady-state 
values, reserving the unprimed status for the quantities referring 
to the •'base” conical flow. To lust order in small quantities, 
equation (II) becomes 

D ! DID 

n, . -- V*//' + Vll'-X 111 P ’(IS) 

l>l n - Dl 

because, in the steady-state conical flow, the stagnation enthldjiy 
is everywhere uniform. This last equation may be further simpli- 
fied by noting that 

In (pir 2 ) — const (10) 

Then, for spherically s\ lumetric flow 


£ ■= inoditied radial variable 


■I 


M 




ir = dimensionless perturbed pressure, p' ip 
rr = also, 11.1 1 15!) 
p -- fluid density 
t dimensionless lime, uA/U 

if) -- phase angle 
a’ •== angular velocity of oseillat. inn 
If solid angle subtended by cone 

Noik: l’rimes denote perturbed variables. 
Subscript zem refers to conditions far upstream. 
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v 2 //' + V//'-V 111 P 


I 0 >•'■!! ' dir id In I- 2 \ 

i 2 d r 0/ dr I dr r f 
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with 
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Tilt partial diffcrel.ind equation fur the Ihictuation of t In* stagna- 
tion enthalpy is. lii- ri fori', 


/*’ ( - I exp | iir i )(£) (2d) 
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( IS) 'I’ln* corresponding differential equation for »? in 
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II is convenient lo romliT equation (10) dinioiisionless by in- 
trodneing Mir follow ing new variables: 
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(21) 
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Acre C: i* (hi* seoustie velocity for tin* tloiii in the reservoir end A* 
is the rmi'.ns ,m< asitrod from I ho apex of I ho coin*) w horo a Mach 
number of unity would bo ohtainoil if (ho conical (low- actually 
persisted up to that railius. In terms oi those variables, equal ion 
(18) becomes 


liquation (till is in a form to which the well-known V KTt 
method |2| can bo applied. The solutions so derived are, in this 
ease, asymptotically exact, in the limits of small Mach number, 
M, or of liiyli win o number, ir. In fact, as shown in the Appendix, 
the error of the \\ Kit method as used here is OfM*/ 12 '*)- The ap- 
proximate solution of equation ttSI ) uliioli has a form at C “■ o> 

eorreBpowliny lo an outwardly radiation wave is 


/> 

/ir 


/ tie \ 2 DX ( r \ o /u„\ d.n' 
\ a / Dr “ df V r / dj- 
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Further progress is aided In performing a, freipiencv analysis in 
time. Therefore we assume that .V can be represented in (lie 
form: 


?; v v * 1 

With Um NiuiK' order of iireui a<\v \\t‘ him* 

W - ^ j '"-I’ j'Act r ■- £) -I fie j Miff -I- r(£)i/£ j- 


(22) 


■v =- ray 


( 22 ) 


Upon insertion of this form for .1C in equation ( 18), we obtain tlm 
following ordinary differential equation for F: 

r (icdF 

' Li - m 
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Amplitude and Phase Relations Between Flow and 
Pressure Distui bailees 

From equations (7) and (tl) we obtain for the disturbed pres- 
sure and velocity the following expressions: 


w here 
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WKB Solution lor ,v 

We now.proceeil to obtain solutions to equation (22). Several 
transformations of variable are necessary to put this equation in 
more tractable form. First, let 


Moreover the. isentropie. relation between pressure and den>.itv 
gives 


Therefore, the perturbed mass flux is 


r 
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Next, let 
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ln dimensionless form, equation (26) becomes 

1_ d7T _ 0.1C M dX 

y,ft* dr dr ^ 1 — M 2 df 
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Correspondingly, equation Bib) heroines 
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The W KB form of .10, given l>y equal ion (ill), i.< noil adapted to 
evaluating this last expression. After some simplification, the 
following form is found; ilms 
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When expressions ( 15), (-Id) mul (51) are used in eipiatiou (4 I), 
the following lion I formula for the nozzle performance is obtained : 


Or 

Or 

(pr)' 

l 

PL. 

Or 


Since, for complex periodic lime variation of rr and (pr)', their 
absolute ratio and their phase are the same as those, of their time 
derivatives, the right-hand side of epilation (4 1) in polar form 
provides the desired answer. We have yet, however, to express 
the answer in terms of readily interpreted physical quantities. 

Formulation in Terms of Strouhal and Mach Numbers 

The following two results can be established alter considerable 
tedious algebraic manipulation: 


nt. 


7,M‘* 

1 ~ M 3 


M 3 + 


i S* 

-f *M -js* 


1 (M)f 

~S* f 




Here ,t(M) is an abbreviation for 


-KM) * 


57* 


3 ,■ (7* - D C*7* - 5) 


M 3 


(54) 


Discussion of Results 

Although entirely accurate only in the limits of low Mach 
number mid high Strouhal number, equation (511) does, never- 
theless, yield the rorieet. quasi-static relations as S* — *• 0; i:e., 
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We observe, among other things, that the upstream influence of 
pressure Hurt nations vanishes as ihe throat Mach number ap- 
proaches unity. 

Now if we assume that •• . ■ 


Those results are valid for the arbiti.irv smg!e-|ihase fluid, exeejit 
that in obtaining equation ( UD the derivative of 7* at constant 
entropy has been neglected. For the monatomic gas, i(£) reduces 
to an especially simple oxpiossion. 

We now proceed to eliminate £f) in the foregoing expressions 
in favor of the Strouhal number. This elimination is accomplished 
ns follows. The Strouhal number at a given point in the flow is 
conventionally defined as ulh'r. To account for the curvature of 
tin- flow area, we here modify tliK definition slightly, taking S -- 
, as 4.4 tr)’ 1 '-! 'J'o relate u'tn i/“, we note that 
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Then, with (l denoting the solid angle of the cone, ivc obtain 
easily enough 


In Figs. 2 and .'!, approximate values of the amplitude factor -V 
mid of the phase angle as found from equation (511), (ire shown 
for various Mach numbers in the ease of a perfect gas with a 
spei ifie heat ratio of 1,4. 

From equal ion (5:1) it, can be seen that the phase angle for M 
■ 0 starts out ;d ISO deg for S + - ~ 0 and approaches .1)0 deg for 
S* — At Hie same time the amplitude factor for M = 0 in- 
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. creases stcuiiilv from unity to infinity. It. is permissible to think 
of the M = 0 condition ns that for an incompressible fluid, for 
■which every portion of t ho thud within the nozzle must oscillate 
in phase. The foregoing phase and amplitude relations reflect, 
then, the increasing difficulty for high-frequency pressure oscilla- 
tions to overcome the inertia of all. the fluid in fliis synchronous 
column within the nozzle. 

When the Mach number at the nozzle is finite, equation (fill) 
shows that 

S* 0 (fill) 

AV’*' -*~ld : S* — a. ((10) 

M 

‘Again, nv obtain quasi-static behavior at low Strouhal numhers. 
The result fof S* -* •» can he expressed in complete form as 


jh<; . 

Or ir 

Xfir)'~ 1 i- M ' 
Or 


(hi) 


A special case of this last formula is obtained when M is small. 
In this ease, the pressure fluet nat ion and mass-llux lluet nation aro 
mutually related in the same maimer ad they are when a plane 
v ail vibrates adjacent to a semi-infinite medium having a finite 
sound velocity. The similarity comes about because, at, high fre- 
quency, the waves within the nozzle are so crowded together that 
the effects of curvature of the phase fronts arc negligible. 

The trends just discussed arc well illustrated by the curves in 
Figs, 2 and 3. In addition, we observe that,, for Strouhal numhers 
If-MS than 0.1, deviations in phase and amplitude from quasi-static 
behavior arc small, and, for Strouhal numbers less than 1.0, de- 
viations between incompressible (M (I) and compressible 
(M > 0) behavior are small. 


One of the conclusions of the present analysis is that the angle 
of nozzle divergence is of some importance. However, for nozzles 
of reasonable design, the range of total included angle is some- 
what, limited, a total included angle of GO deg being representa- 
tive. For this particular angle, the modified and usual Strouhal 
numbers are identical. , 
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APPENDIX 

The purpose of this Appendix is to ascertain the order of mag- 
.. nifudeof the error incurred through use of the lV r KIl method iri 
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the Solution of equation (;il). The treatment employed will 
follow completely, for our specific case, the analysis given in 
Section 3.3 of reference |4). 

We have the differential equation 
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- r(4)jt)(4) - 0 

«•>•->) 
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Now let. 



V(0 ~ 

= c -’"f:(4) 

(tin) 

The resulting dilTerential equation for 2 is, then 


2 "(4) “ 2fm? 

'(4) - '<£> - « 

((Mi) 

A solution of the foregoing difterential equation is provided hy 
that 2 ( 4 ) which satisfies the following integral equation: 

2(4) - 1+1' J 

— A /.(|) : (,)rvi 

2ie 

(f>7) 

««)*= -4V(4); |/.|</> 

(t>S) 


2 ( 4 ) in equation (07) can he expressed by the infinite series: 


and it. can he seen that the series (0!)) converges for all finite 
£ and tc. , : 

As a special rn.se of (70), v e have 

Ulylt -1- -■ f r(l)<lt (72) 

2»’J f 



Now r(f ) is zero at. infinity, and we shall presume that it increases 
inonotonically with decreasing Then, through integration hy 
parts, we can show that 
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As 4 approaches infinity, it is evident that 
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Then 


r ™ OtM- 1 ); 1) = 0(M\) 
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(70) 1 


Reference now to equation (34) of the text will show that the 
WKli approximation for 2(4) is '* 

■ f l r(t)tlf 

2(4) Sf'' • ; ■ 


*($> » E 

0 


where 2 , 4 4 ) is defined to he unity, and 
j | _ 


2»+i($) = ' j 


Ul)zAt)l~W 


Asa consequence of equation (<i7), it can he proven that 

I />" 


k(4)| < 


«! (i/>4)" 



we find hy int.ereomparison of series (7(i) and (77) that the frac- 


( 71 ) lional error of the WKl! method as employed in this work is 
OCM’/ie*). 


Printed hi II. S. A. 


8 


Transactions of the ASME 



